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Abstract 

We analyze the quantum entanglement between two interacting àtoms trapped in a spherical 
harmònic potential. At ultra-cold temperature, ground state entanglement is generated by the 
dominated s-wave interaction. Based on a regularized pseudo-potential Hamiltonian, we examine 
the quantum entanglement by performing the Schmidt decomposition of low-energy eigenfunctions. 
We indicate how the àtoms are paired and quantify the entanglement as a function of a modified 
s-wave scattering length inside the trap. 
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I. INTRODUCTION 



The interactions between trapped ultra-cold àtoms govern manv interesting collective 
quantum phenomena, ranging from Bose-Einstein condensation l| to recently observed 
fermion superfluid Q|. At sumciently low energies, it is known that short ranged atom-atom 
interactions can be replaced by a pointlike regularized pseudo-potential under the shape 
independent approximation 3]. The strength of such a pseudo-potential is determined by a 
s-wave scattering length a only, and hence one can control atom-atom interaction by tuning 
the scattering length via the technique of Feshbach resonance. For trapped svstems, the 
theory of pseudo-potential has been examined in details by several authors 0,0]. As long 
as the range of the actual atom-atom interaction is short compared with the width of the 
trap, low energy eigenfunctions can be accurately captured by the eigenfunctions of the 
pseudo-potential, except for a few tightly bound states that may exist inside the range of 
the inter-atomic potential. Such a universal applicability is the essence of shape independent 
approximation. Therefore the study of eigenfunctions of pseudo-potentials would provide 
useful insight about genèric features of two-body correlations in the low energy regime. 

In this paper, we address a fundamental question of how the scattering length controls 
quantum correlations between two ultra-cold àtoms inside a harmònic trap. Quantum con- 
trol of trapped ultracold àtoms has been a subject of considerable research interest, regarding 
potential applications in quantum information processing jf], Q|. For example, col·lisions of 
ato,„ 8 can be exp.oited to perfonn vaxious ouantum lògic operations Q. However, the na- 
ture of quantum entanglement arising from s-wave scattering has not been fully explored |8J . 
Such an entanglement is inherent in the continuous degree of freedom of àtoms, and it may 
have effects on the fidelity of quantum gates based upon collisional mechanisms In this 
paper we will analyze the quantum entanglement of the low energy eigenstates defined by 
the regularized pseudo-potential and the harmònic trap. By performing the Schmidt decom- 
position of low energy eigenfunctions, we show that quantum entanglement is manifested 
as pairing of àtoms in a set of orthogonal mode functions in three-dimensional space. In 
particular, the angular momenta are identified as good quantum numbers to characterize 
the Schmidt mode functions. We will present numerical results that quantify the degree of 
entanglement as a function of the scattering length a. In addition, we will examine the en- 
tanglement in the a — > oo limit. Such a strong coupling regime corresponds to the unitarity 
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limit in degenerate quantum gases |9j. The study of pairing in two-body models in such a 
limit may shed light on quantum correlations in the more difRcult many body problems. 



II. REGULARIZED HAMILTONIAN AND ENERGY EIGENSTATES 

To begin with, we consider a system of two interacting àtoms with equal mass trapped 
in a spherical harmònic potential. The Hamiltonian of the system is given by, 

H = -— V? - — Vo + -mu 2 r 2 + -mu 2 r 2 + VYri - r 2 ), (1) 
2m 1 2m 2 2 1 2 2 V ; ' K J 

where ri and r 2 are the position vectors of the two àtoms, m is the mass of the atom, 
and uj is the trap frequency. The interaction between the two àtoms is described by the 
short range potential V which will be replaced by a pseudo-potential in Eq. (3) under the 
shape independent approximation. For convenience, we separate the Hamiltonian into a 
center-of-mass part and a relative part, H = H cm + H re i so that 

H cm = --V]j + 2R 2 , (2) 

o 

Hrei = -W r + \r 2 + 2na5^ (r) |r, (3) 
2 2 or 

with R = (ri + r 2 )/2 and r = r x — r 2 . Here the energy and length are expressed in units of 
Tiio and (h/iiuj) 1 ^ 2 (where fj, — m/2 is the reduced mass) respectively. The strength of the 
pseudo-potential is characterized by the modified s-wave scattering length a. For a given 
inter-atomic potential, the precise value of a depends on the trappjng potential and it can 
be determined self-consistently by the methods discussed in Ref. pj. In this paper we will 
treat parameter of the Hamiltonian. 

The s-wave eigen-functions of the H have been solved analytically in Ref. Q. Given a 
scattering length, the eigen-energy E of H re i is defined by the solution of: 

r(-E/2-l/4) 

2r(-£/2 + 3/4) v ; 

with T being the gamma function. The eigenfunctions of H re i with the energy E are given 
by: 

^(r) = Ae~ r2 / 2 U (-E/2 + 3/4, 3/2, r 2 ) (5) 



where A is a normalization constant, and U(a,j3,z) is the Kummer's function 12|. 



In this paper we assume that the center of mass wave function is the ground state of 
H cm , which is a simple gaussian: $(R) = 2v / 2e~ 2i?2 /V 3 / 4 . Combining $(R) with ^(r), the 
two-particle energy eigenfunctions are given by 

^(r 1; r 2 ) = $(R)^(r). (6) 

In Fig. 1, we illustrate how the eigen-energies depend on the scattering length. For conve- 
nience, we choose to plot with the inverse of scattering length, i.e., 1/a, in order to indicate 
the continuous branch associated with the ground states. It should be noted that the ground 
state energy becomes large and negative when 1/a is positive and large. This feature also 
occurs in the absence of the harmònic trapping potential, and it is due to the existence of a 
tightly bound state in the a — > + limit As an illustration, we show in Fig. 2 the radial 
probability density associated with relative coordinate wave functions at several vàlues of a. 



III. SCHMIDT DECOMPOSITION 



The characterization of quantum entanglement is achieved by Schmidt decomposition of 
^ (i - !, r 2 ), which reads: 

*(ri,r2)=E^(ri)«i(r2), (7) 
i 

where \j are eigenvalues, and Uj and Vj are Schmidt eigenmodes defined by, 

J dr[ J dr 2 Ï! (n, r 2 ) ^* (r' 1; r 2 ) Uj (r[) = À^n), (8) 
J dr' 2 J dv^ (r 1; r 2 ) V* (n, r' 2 ) Vj (r' 2 ) = X jVj (r 2 ). (9) 

Note that the mode functions Uj form a complete and orthonormal set, and the same is true 
for Vj. If the atom 1 appears in the mode Uj, then with certainty the atom 2 must be in the 
mode Vj. In other words, Eq. (7) indicates the pairing structure of the two particle state. 
In addition, the distribution of provi des a measure of the degree of entanglement. This 
is usually discussed in terms of the entanglement entropy S = —J2j^j^og\j. However, a 
more transparent measure is the effective number of Schmidt modes, which is provided by 
the Schmidt number: K = l/J^jX'j 10]. A disentangled (product) state corresponds to 
K = 1, i.e., there is only one term in the Schmidt decomposition. The larger the value of 
K, the higher the entanglement. We point out that 1/K equals the purity of the density 



matrix of an individual particle. The purity has also been employed as a measure of the 



degree of entanglement in various physical situations 



ui 



To carry out the Schmidt decomposition of the wave functions in three dimension, we 



note that R = yrj + r 2 + 2r x r 2 cos 7/2 and r = y r\ + r\ — 2r x r 2 cos 7, where 7 is the 

angle between 17 and r 2 . Therefore the wave function \l/(r 1 ,r 2 ) = \& (77, r 2 , cos 7) = 

00 

on (77, r 2 ) Pi (cos 7), where Pi(x) is the Legendre polynominal, and 



1 r n 

— / G?7sin7 ■ \l/ (77, r 2 , 00S7) Pi (cos 7). (10) 
Jo 



2/ 

(ri,r 2 ) = — 

z JO 

With the help of the addition formula: (2/ + l)P t (cos 7) = 4rr £ (0 1} X ) Y" ím (0 2 , 2 ), 

m=— i 

we have 

*(r 1) r 2 ) = 47rf; É ^^y^^.^y^^,^). (11) 

1=0 m=-l 

This expression is already partially in the Schmidt form, because the pairing of angular 
functions has been identified. The remaining task is to decompose ai (77, r 2 ) for each l. This 
can be achieved by performing the Schmidt decomposition of the function rir 2 «/(ri, r 2 ), i.e., 



r\r 2 OLi (ri, r 2 ) = J2\ ^mu n i (n) v ni (r 2 ). (12) 

n 

Here the prefactor r^r^ is introduced in order to ensure correct normalization in radial 
directions. 

The final form of the Schmidt decomposition of the wave function Eq. (5) now reads: 

00 00 1 



*(ri,r 2 ) = ££ E 

n=l 1=0 m=—l 



21 + 1 77 



r 2 



(13) 



Our derivation indicates a general feature that for any wave function that is a function 
of distances R and r only, the angular parts of Schmidt modes are simply the spherical 
harmonies. The pairing involves angular momentum quantum numbers (l,m) and (l,—m), 
i.e., the same / and opposite m must be paired. In addition, each m is of equal weight for 
a given /. Therefore if one could select Schmidt modes with a fixed n and / via projective 
measurement, then the resulting state is a maximally entangled state among various m's on 
the projected / manifold. 
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IV. NUMERICAL RESULTS 



The fact that the angular parts of Schmidt modes can be obtained analytically reduces 
the còmput ational difficulty in the original three-dimensional problem. Finding u n i, v n i and 
X n i in Eq. (12) is a relatively simple numerical task, because r 1 r 2 ai(ri,r 2 ) behaves as a 
two-particle (one dimensional) wave function in non-negative r% and r 2 regions. Specifically, 
we first obtain a.i(ri,r 2 ) by performing numerical integration of (10) for discretized vàlues 
of ri and r 2 , typically with the spacing Ar = 0.01. For the s-wave eigenfunctions considered 
here, it is sufficient to choose r± and r 2 ranging from to 3.5, where the wave functions 
are mainly confined. Because of the symmetry properties of r 1 r 2 a;(r 1 , r 2 ), u n \ and v n \ are 
the same real functions. Therefore u n i, v n \ and X n i can be obtained by diagonalizing the 
matrix rir 2 ai(ri, r 2 ). For low energy states considered in this paper, l up to 30 are typically 
sufficient in order to obtain convergent results. 

We can now obtain the value of K for the state (13) from the Schmidt eigenvalues, i.e., 

oo oo 

^ = VEEA^. (14) 

n=l 1=0 

Here A ni = 16vr 2 À n ,/(2Z + l) 3 / 2 is defined. The main result of this paper is shown in Fig. 
3, where the vàlues of K are displayed as a function of the inverse of the scaled scattering 
length for low energy states. We notice that higher excited states generally have higher 
quantum entanglement. However, the ground state (curve A) shows a distinct behavior in 
the positive 1/a region, where we notice a sharp rise of K as 1/a increases. 

The strong ground state entanglement in the large positive 1/a limit is understood as the 
appearance of increasingly bounded àtoms. This is implied in the ground state energy curve 
(A) in Fig. 1, as well as in the wave function (1/a = 2) shown in Fig. 2. A crude estimation 
of K can be made from the analytical results of gaussian functions. For gaussian functions 
separable in the center of mass and relative coordinates, it is known that K oc (Ai?) 3 / (Ar) 3 
when the center of mass width Ai? is much wider than the width of the relative coordinate 
Ar 14]. Here the tightly bound state corresponds to a strong localization in particle's 
relative distance, with Ar of order a in the 1/a ^> 1 limit. Therefore AR/Ar ^> 1, and 
hence high vàlues of K can be expected as in the case of gaussians functions. However, we 
remark that gaussian can only serve as a guide here, because the singular 1/r dependence 
in ^e{t) cannot be captured by gaussians. Indeed, K depends on a in a complicated form 
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according to our numerical calculations. 

For excited eigenstates on the branches B and C, we see an interesting feature in Fig. 3 
that the change of quantum entanglement is only sensitive to a range of scaled scattering 
lengths. Such a window of scaled scattering lengths is highlighted by dash lines in Fig. 
3, where we found that K changes significantly with 1/a when l/\a\ < 2. Recalling that 
we are using the length unit defined by the harmònic trap, our results suggest that the 
s-wave interaction can appreciably affect quantum entanglement in excited states when a is 
greater than or comparable with the width of a ground state particle inside the trap, i.e., 

1 /2 

\a\ > (h/2mw) 1 . This feature is also true for the ground state with negative scattering 



The structure of quantum entanglement is characterized by the distribution of Schmidt 
eigenvalues and the Schmidt mode functions. In Fig. 4 we show the distribution A n i for 
the ground state with 1/a = —2 and 1/a = 2. In the case of 1/a = —2, the entanglement 
K — 1.2 is not high because of the presence of a dominant (n, l,m) = (1, 0, 0) Schmidt mode, 
which covers about 70% probability of the state. On the contrary, K = 10.13 is much higher 
for the 1/a = 2 case. This is shown in the distribution of A n i (Fig. 5b) in which Schmidt 
modes with n = 1 and higher angular momentum number / are involved. In other words, 
the strong entanglement is mainly manifested in the angular variables. 

In Fig. 5, we illustrate the shape of leading Schmidt modes corresponding to the ground 
state wave functions in Fig. 4. Since the angular part are simple spherical harmònic func- 
tions, we show only the radial mode functions in the figure. Apart from the fact that 
Schmidt modes of positive scattering lengths are more attracted to the origin, the change of 
scattering length has small effects on the mode functions. This is in contrast to the Schmidt 
eigenvalues, which are more sensitive to the vàlues of a as depicted in Fig. 4. 

Finally, let us discuss the limit of jaj — > oo or l/|a| — > 0. In such a limit, we find that the 
eigenfunctions take simple analytic forms. We present some of the s-wave eigenfunctions of 
H in the limit lol — > oo: 



lengths. 



(ri,r 2 ) 



2 e" r ?e- r 2 



(15) 



*n (ri,r 2 ) 



71-3/2 | Fl _ r2 | ' 

y/2 e- r h~ r ï - 

n 3/2 | ri _ r2 | . 



1 - (ri - r 2 f , 



(16) 



^12 (ri,r 2 ) 



e-^e-i r 



l-4(r 1 -r 2 ) 2 + ^(r 1 -r 2 ) 4 . 



(17) 



7r 3 / 2 |ri — r 2 
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Here Eqs. (15-17) correspond to the energy eigenfunctions (center of mass + relative coordi- 
nates) associated with the lowest three states of H re i, and the ground state of H cm . Although 
Eqs. (15-17) are simple expressions, the Schmidt decomposition still cannot be carried out 
analytically. We perform numerical calculations which give: K w = 1.98, K\\ = 3.45, and 
K12 = 9.11 for the wave functions given in Eqs. (15-17). Therefore the degree of entangle- 
ment remains finite at the infmite scattering limit. 

We remark that the regularized pseudo-potential method can only describe wave functions 
at inter-atomic distance much larger than the range of the actual inter-atomic potential b, 
i.e., |ri — ^> b. Inside the interaction range b, the actual wave functions remain finite 
as |ri — r%\ — > 0. Therefore the singular behavior at ri = r 2 in Eqs. (15-17) is only an 
artifact of the shape independent approximation, and hence these wave functions should be 

1 /2 

understood for |ri — r2 1 b only. However, since b is typically much smaller than (K/fiu)) 
(which is about the width of 1 í r ), the probability of finding the two partides within the range 
b is negligible Q|. This justifies the use of the shape independent approximation here. 



V. SUMMARY 



To summarize, we present a procedure to analyze the s-wave quantum entanglement 
between two ultra-cold àtoms in a spherical harmònic trap. The s-wave interaction is de- 
scribed by the regularized pseudo-potential. By performing the Schmidt decomposition of 
low energy eigenstates, we quantify the quantum entanglement and indicate its dependence 
on the modified scattering length. In particular, our Schmidt analysis reveals the angular 
correlations by showing explicitly the pairing of spherical harmònic functions. For small and 
positive a, the ground states are highly entangled states, and we explain this feature as a 
consequence of tight binding between the àtoms. For low excited states and ground states 
with negative a, we find that the atom-atom interaction can only appreciably affect the 
entanglement when the scattering length is larger than the width of the (non-interacting) 
ground state defined by the trap. However, the degree of entanglement remains finite in 
the large scattering length limit. Our work here indicates that quantum entanglement can 
be controlled by the scattering length. To explore applications of s-wave entanglement, one 
may need to establish schemes for detecting Schmidt modes. In addition, the dynamics of 
entanglement associated with non-stationary states of the system is also an interesting tòpic 
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for open future investigations. 
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FIG. 1: Eigen-energies of H re i as a function of the inverse of dimensionless scaled scattering length. 
E is in units of fico. The three branches: A, B, and C correspond to the lowest three states with 
zero angular momentum (in relative coordinate). 
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FIG. 2: Radial densities p(r) = Airr 2 \ipE{r)\ 2 associated with the ground state wave functions at 
l/a = 0, ±2. The radial distance r is in units of (h/ ^lo) 1 ^ 2 . 
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FIG. 3: Schmidt number as a function of the inverse of dimensionless scaled scattering length for 
the low eigenstates associated with the three energy curves in Fig. 1. K is plotted up to 1/a = 2 
for the ground state curve A because the high vàlues of K is out of the range of the figure for 
1/a > 2. 




FIG. 4: (Color online) Distribution of A n ; (see text) of ground states for: (a) 1/a = —2 and (b) 
1/a = 2. 
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FIG. 5: Probability density of the radial part of Schmidt modes p n \ = |u n ;| 2 , where the u n \ is 
defined in Eq. (12). The sòlid line corresponds to the ground state with 1/a = —2, while the dash 
line corresponds to 1/a = 2. The radial distance r is in units of (h/ fiuj) 1 ^ 2 . 
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